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Artificial molecules containing just one or two electrons provide a powerful platform for 
studies of orbital and spin quantum dynamics in nanoscale devices. A well-known example of 
these dynamics is tunnelling of electrons between two coupled quantum dots triggered by 
microwave irradiation. So far, these tunnelling processes have been treated as electric-dipole- 
allowed spin-conserving events. Here we report that microwaves can also excite tunnelling 
transitions between states with different spin. We show that the dominant mechanism 
responsible for violation of spin conservation is the spin-orbit interaction. These transitions 
make it possible to perform detailed microwave spectroscopy of the molecular spin states of 
an artificial hydrogen molecule and open up the possibility of realizing full quantum control 
of a two-spin system through microwave excitation. 
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n recent years, artificial molecules in mesoscopic systems have 
drawn much attention owing to a fundamental interest in their 
quantum properties and their potential for quantum informa- 
tion applications. Arguably, the most flexible and tunable artificial 
molecule consists of coupled semiconductor quantum dots that are 
defined in a two-dimensional electron gas using a set of patterned 
electrostatic depletion gates. Electron spins in such quantum dots 
exhibit coherence times up to 200 jis (ref. 1), about 10^-10^ times 
longer than the relevant quantum gate operations'"^, making them 
attractive quantum bit (qubit) systems.^ 

The molecular orbital structure of these artificial quantum 
objects can be probed spectroscopically by microwave modulation 
of the voltage applied to one of the gates that define the dots^. In 
this way, the delocalized nature of the electronic eigenstates of an 
artificial hydrogen-like molecule was observed^'^. More recently, 
electrical microwave excitation was used for spectroscopy of single 
spins^"^^ and coherent single -spin controP'^^ through electric- dip ole 
spin resonance. 

Here we perform microwave spectroscopy^'^'^' on molecular spin 
states in an artificial hydrogen molecule formed by a double quan- 
tum dot (DD) which contains exactly two electrons. In contrast to all 
previous photon-assisted tunnelling (PAT) experiments, we observe 
not only the usual spin- conserving tunnel transitions, but also tran- 
sitions between molecular states with different spin quantum num- 
bers. We discuss several possible mechanisms and conclude from 
our analysis that these transitions become allowed predominantly 
through spin-orbit (SO) interaction. The possibility to excite spin- 
flip tunnelling transitions lifts existing restrictions in our thinking 
about quantum control and detection of spins in quantum dots, and 
allows universal control of spin qubits without gate-voltage pulses. 

Results 

Device and excitation protocol. Figure la displays a scanning 
electron micrograph of a sample similar to that used in the 
experiments. It shows the metal gate pattern that electrostatically 
defines a DD and a quantum point contact (QPC) within a GaAs/ 
(Al,Ga)As two-dimensional electron gas. An on-chip Co micro- 
magnet (|a magnet) indicated in blue in Figure la generates an 
inhomogeneous magnetic field across the DD, which adds to the 
homogeneous external in-plane magnetic field B (Methods), but is 
not needed for the molecular spin spectroscopy. The sample was 
mounted in a dilution refrigerator equipped with high-frequency 
lines. The gate voltages are set so that the DD can be considered 
as a closed system (the interdot tunnelling rates are 10^ times 
larger than the dot-to-lead tunnelling rates), and the tilt of the 
DD potential is tuned by the dc-voltages and V^, applied to 
the left and right side gates. Working near the turn- on of the first 
conductance plateau, the current through the QPC, /gp^, depends 
upon the local charge configuration and provides a sensitive meter 
for the absolute number of electrons {n^, n^) in the left and right 
dot, respectively^^'^^. 

First, we excite the DD as indicated in the top panel of Figure 
lb, by adding to continuous-wave microwave excitation at fixed 
frequency v= 11 GHz. When the photon energy of the microwaves 
matches the energy splitting between the ground state and a state 
with a different charge configuration, a new steady-state charge 
configuration results, which is visible as a change in the QPC cur- 
rent, AJqpc. The excitation is on-off modulated at 880 Hz and lock- 
in detection of A/qp^ reveals the microwave -induced change of the 
charge configuration (see Methods for further experimental details). 
The lower panel of Figure lb shows A/gp^ as a function of and 

near the (1,1) to (0,2) boundary of the charge stability diagram. 
Sharp red (blue) lines indicate microwave -induced tunnelling of 
an electron from the right to the left dot (left to right), labelled as 
An = + 1 (An = - 1), respectively (Fig. Ic). Sidebands can result from 
multi-photon absorption. At the boundaries with the (0,1) and (1,2) 
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Figure 1 1 Photon-assisted tunnelling in a 2-electron double quantum 
dot. (a) Scanning-electron micrograph top view of the double-dot gate 
structure with Co micromagnet (blue). The voltages applied to the left 
Vl and right side gates (red) control the detuning e of the double-dot 
potential. The double-dot charge state is read out by means of the current 
/qpc running through a nearby quantum point contact (white arrow), 
(b) Charge stability diagram around the 2-electron regime at 6 = 1.5 T. 
(riyn^) indicate the absolute numbers of electrons in the left and right 
dot, respectively. During measurements, 11 GHz microwaves with 880 Hz 
on-off modulation are applied to the right side gate. The top panel displays 
schematically one cycle of the ac signal applied to the right side gate. 
Along the detuning axis (dashed arrow) PAT-lines are observed, (c) In the 
conventional picture of PAT, the first sidebands seen in b should appear 
when the detuning of the (0,2) and (1,1) states matches the photon energy, 
and interdot tunnelling is induced. Further sidebands are then interpreted 
as multi-photon transitions, jip is the chemical potential of the left and right 
electron reservoir, (d) Same as in b, but the microwaves are interrupted 
every 5jis by a 200ns, P^ = 2m\/ detuning pulse applied to the left and 
right side gates (see the schematic in the top panel). The pulses generate 
a reference line (black arrow) due to mixing at the ST^ anti-crossing. 



charge states, no energy quantization is observed, because, here, 
electrons tunnel to and from the electron- state continuum of the 
leads. At first sight, the observations in Figure lb thus appear to be 
well explained by the usual spin- conserving PAT processes. 

Surprisingly, the position in gate voltage of the resonant lines 
exhibits a striking dependence on the in-plane magnetic field, 
B. This is clearly seen in Figure 2a,b, which display the measured 
PAT spectrum along the DD detuning e axis (dashed black arrow in 
Fig. lb) as a function of B for 20 and 1 1 GHz excitation, respectively. 

As the gate constitutes an open-ended termination of the trans- 
mission line, the excitation produces negligible AC magnetic fields 
at the DD (estimated in Methods), and is, therefore, expected to give 
rise to only electric -dip ole -allowed spin-conserving transitions, 
with no B dependence. Furthermore, there is a pronounced asym- 
metry between the position of the red and blue PAT lines. 

In these figures, the detuning axis was calibrated for all mag- 
netic fields by introducing a reference line (Fig. Id, black arrow) 
that facilitates interpretation of the spectra despite residual orbital 
effects of the magnetic field. This line was produced by interspersing 
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Figure 2 | Photon-assisted-tunnelling spectra and simulations. (a,b) Microwave-induced change of the QPC current A/gp^ as a function of the double- 
dot detuning e and the external magnetic field 6 for 20 GHz (a) and 11 GHz (b) frequency, respectively. Singlet-triplet mixing due to 2-mV detuning pulses 
generates a reference signal that is used to calibrate the detuning axis (see lower panel in Fig. 2f). (c) Eigenenergies versus double-dot detuning e of the 
two-electron spin states in the (1,1) and (0,1) charge regime for two external magnetic fields 6 = 2. 5T (upper panel) and 6 = 1. 5T (lower panel), respectively. 
5(0,2), 5(1,1) and 7(1,1) character of the eigenstates is indicated by blue, green and red colour, respectively. The molecular spin ground state is indicated by 
thick lines. The vertical arrows indicate PAT transitions for a constant microwave frequency involving spin flips. The transition indicated by a dashed arrow 
is suppressed, because the initial state lies above the ground state. The red circle in the lower panel indicates the detuning position of the reference signal, 
that is generated by a detuning pulse with amplitude to the 57^ anti-crossing. (d,e) Simulated PAT spectra for 20 GHz (a) and 11 GHz (b) frequency, 
respectively. The colour indicates the change of the population of the steady-state charge state An as would be observed in an on-off lock-in detection. A 
finite temperature of 100 mK and spontaneous relaxation through the phonon bath are taken into account, (f) /qp^- scanned with higher resolution in the 
anti-crossing region (black rectangle in Fig. 2b) at 11 GHz excitation. The green dashed lines indicate the expected detuning positions of the PAT transitions. 
The horizontal black dashed line indicates the magnetic field, at which the electron spin resonance condition is fulfilled E^-giiQ{?> + bQ)-hv. The graph is 
concatenated from two scans that overlap at 1.9 T. The inset displays three magnetic fields (blue dots), at which the centre of the horizontal blue triplet 
resonance line is observed, as a function of the microwave frequency v. The dashed black line gives the expected position of the triplet resonance. 



the microwaves every 5 \is with 200 ns gate-voltage pulses along the 
detuning axis (see the top panel in Fig. Id), leading to singlet-triplet 
mixing as described in ref. 2. The short gate-voltage pulses do not 
noticeably alter the position of the PAT lines (compare Fig. lb,d). 
The reference peaks visible at around e=200|aeV in Figure 2a,b 
were aligned by shifting all data points at a given B by the same 
amount in detuning (see Supplementary Note 1 'Calibration of the 
detuning axis' for the full details of this post-processing step and 
Supplementary Fig. Sla-d for the corresponding spectra). 

Interpretation of the photon-assisted tunnelling spectra. The 

complexity of the PAT spectra shown in Figure 2a,b can be under- 
stood in detail, if we allow for non-spin conserving transitions. 
The two diagrams in Figure 2c show the energies of all relevant 
DD-states (four (1,1) -states and one singlet S(0,2) -state) as a func- 
tion of £ for two different (fixed) magnetic fields, that is, the spec- 
trum of the DD along the two horizontal dotted lines in Figure 2b 
Note that the only difference between the two diagrams is the split- 
ting between the three triplet T( 1,1) -states. 

First, we explain the resonances observed along the upper dotted 
line in Figure 2b (5 = 2.5T). In the corresponding (upper) diagram 
in Figure 2c, we plot the ground state energy for all e with a thick 
line. If the microwave excitation is off-resonance with all transitions, 
the system will be in this ground state. For instance, at 150|ieV, 
there is no state available 1 1 GHz above the S(0,2) ground state (grey 



arrow) and the system stays in S(0,2). However, when decreasing the 
detuning, at some point T+(l,l) becomes energetically accessible 
(red arrow) and, because we allow for non- spin- conserving transi- 
tions, is populated owing to the microwave excitation. For this PAT 
transition, the spin projection on the quantization axis is changed 
by Am = -i- 1. The resulting change of steady-state charge population 
(increased population of (1,1), or A/qpc>0) is detected by the QPC 
and yields the red peak in Figure 2b. Decreasing e further, there 
are two more resonances detectable: (i) the S(0,2) - S(l,l) transition 
(dotted red arrow, Am = 0), although the signal will be weakened 
owing to the fact that S(0,2) is not unambiguously the ground state 
anymore. Note that a transition S(0,2) - T^(l,l) could appear in 
nearly the same detuning position, as will be discussed below, (ii) 
the T+(l,l)-S transition (blue arrow), where S stands for the hybrid- 
ized S(0,2) -S(l,l) singlet, results in a negative (blue, A/qpc<0, 
Am = - 1) signal from the charge detector as the ground state is now 
(1,1) and the excited state is partly (0,2). We see that this simple 
analysis explains both the positions and the signs of the resonances 
observed in the data. 

A similar analysis can be made for other magnetic fields. For 
instance, for the spectrum plotted in the lower diagram of Figure 
2c, we find two resonances with A/qpc>0 (red arrows), and one with 
A/qpc<0 (blue arrow). Note that the 'blue' transition now connects 
the ground state to the other branch of the hybridized S compared 
to the high magnetic field case. Indeed, the singlet anti-crossing is 
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Figure 3 | Analysis of the photon-assisted-tunnelling spectra, (a) The 

detuning difference Ae between the spin-conserving line Am = 0 and the 
Am = ±1 line is plotted as a function of the external magnetic field B at 
20 GHz (see Fig. 2a), to fit the absolute effective electron g-factor g from 
the slopes of the linear fits (solid lines), inset, Aefor the Am= +1 line at 
11 GHz for different tunnel couplings. The offset of the curves increases 
with increasing (red to black circles) while g remains constant, (b) Ae 
between the Am = 0 line and the Am= -1 (1,1) to (0,2) transition from 
Figure 2c. The fit of the anti-crossing (red line) allows for a precise 
determination of t,. 



directly probed, resulting in the two blue curved lines observed in 
the data around £ = 0 (Fig. 2b). The fading out of the blue signal 
at low fields can be understood from pumping into the metast- 
able state S( 1,1): the microwaves excite the system from T+(l,l) to 
S(0,2), from where it relaxes quickly to S(l,l). However, relaxation 
from S(l,l) back to the ground state is slow because of the small 
energy difference of this transition and the small phonon density of 
states at low energies^^'^^. This pumping weakens the detector signal, 
as S(l,l) has the same charge configuration as the ground state. 

To verify this interpretation, we calculate in Figure 2d,e the 
position and intensity of the spectral lines at fixed microwave 
frequency, on the basis of the energy level diagram of Figure 2c. In 
the simulations, all single-photon transitions between the ground 
state and the excited states are allowed by including a matrix ele- 
ment |T-(1,1)XS(0,2)| (Methods). The input parameters for the cal- 
culation of the resonant positions are the interdot S(l,l) to S(0,2) 
tunnel coupling t^, the absolute electron ^-factor |^|, a magnetic-field 
contribution from the |l magnet parallel to B as well as a magnetic 
field gradient AB = (AB^ ,0,AB^ between the dots (in Fig. 3, we 
show how t,, 1^1 and Bq can be extracted from the experimental spec- 
tra). The colour scale represents the calculated steady- state An that 
results from microwave excitation, orbital hybridization, and pho- 
non absorption and emission at 100 mK. All the PAT transitions vis- 
ible in the simulation also appear in the experiment, with excellent 
agreement in both the position and relative intensity of the spectral 
lines. Especially, the vanishing signal due to spin pumping is also 
predicted by the simulations that include phonon relaxation. Sur- 
prising at first glance, the intensities of the spin-flip lines Am =±1 
and the spin- conserved line Am = 0 are about the same, which holds 
true for both simulation and measurement, although Am = 0 has a 
larger matrix element. This is because the same matrix elements 
enter the spontaneous phonon-mediated relaxation rates, resulting 
in similar steady-state populations for the two transitions. 

When we zoom in on the boxed region of Figure 2b, we see an 
extra horizontal blue feature at 5-2 T (Fig. 2f ) that also appears in 
the calculated spectra of Figure 2e. This feature is due to a triplet 
resonance from T+( 1,1) to 7^(1,1) that becomes detectable by relax- 
ation into the meta-stable S(0,2) state. In the detuning range where 
this line appears, the S(0,2) state lies energetically only slightly above 
the T+(l,l) state, so relaxation back to the T+(l,l) ground state 
is suppressed, again by the small phonon density of states at low ener- 
gies (see Supplementary Fig. S2a for a schematic energy diagram). 
The triplet resonance is expected to appear at E^=giig(B + bQ) = hVy 



where h is Planck's constant and /ig the Bohr magneton. The inset of 
Figure 2f shows the magnetic fields corresponding to the centre of 
the measured triplet resonance line for three excitation frequencies 
(see Supplementary Note 2 'Triplet spin resonance' for details and 
Supplementary Fig. S2b,c for extra spectra taken at different micro- 
wave frequencies), which are in good agreement with the expected 
positions (black dashed lines in Fig. 2f ) based on the values |^| and 
Bq determined in the next section from other features of the PAT 
spectra. Surprisingly, the measured triplet resonance exhibits a finite 
slope in the B(£) spectra. A longitudinal magnetic field gradient AB'I 
gives rise to such a detuning dependence, but the ARU required in 
our simulations to reproduce the observed slope is ABl'^SOmT per 
50 nm, an order of magnitude larger than the gradient we calculate 
for the [I magnet. The magnitude of the slope remains a puzzle. 

Extracting artificial molecule parameters. We now show how |^|, 
and the parameters used for all simulations, can be extracted 
independently from the experimental spin-flip PAT spectra. For this 
analysis, we only use the relative distance Ae between PAT lines at fixed 
magnetic field, to be independent from the calibration of the detuning 
axis by means of the reference line. Figure 3a shows Ae- as a function 
of B using the 20 GHz data. Ae- is defined as the difference in detuning 
between the red Am =±1 and Am = 0 PAT lines as plotted in Figure 2a. 
For a fixed V, both Ae- increase linearly with the Zeeman energy and 
therefore allow fitting of |g|. (Note that for Ae^ the linearity is only 
exact for sufficiently large 5, at which the singlet anti- crossing does 
not affect the detuning position of the S(0,2) to 7^(1,1) transition. 
A least-squares fit to the Ae+ data gives |^| = 0.382 ±0.004 (Fig. 3a). 
From the linear behaviour of Ae^, we also deduce that there is most 
likely negligible dynamic nuclear polarization in the experiment. 

Knowing |^| precisely, we make use of the blue anti- crossing in 
Figure 2b,f to determine (and bo). Figure 3b shows the difference 
in detuning Ae between the blue Am = -1 and the red Am = 0 lines 
in Figure 2f Assuming the Am = 0 line corresponds to the S(0,2)- 
S(l,l) transition (as shown below), then 



Ae'^ 



-{hv-gliLBiB + bo)f 

hv-gllBiB + bo) 



-^(hvf-(2t,)\ (1) 



where the first term is the detuning position of the T+(l,l)-S tran- 
sition and the second the one of the S(0,2) -S(l,l) transition. The 
best fits are obtained with t^ = S.7±0.1 |ieV and = 109 ± 16 mT. This 
value for b^ matches our simulations of the stray field of the [i mag- 
net at the DD location very well (see Methods). 

An important question left open so far is whether the red Am = 0 
line involves predominantly transitions from S(0,2) to S(l,l) or to 
T^(l,l). The transition to S(l,l) does not require a change in the 
(total) spin and is thus expected to be excited more strongly than 
that to T^(l,l). However, relaxation from S(l,l) back to S(0,2) will 
be stronger as well, so it is not obvious what steady-state popula- 
tions will result in either case. Furthermore, given the small energy 
difference between S(l,l) and T°(l,l), the two transitions are not 
resolved in Figure 2. Figure 3a helps to answer this question: The 
observation that Ae^ > Ae" indicates that the Am = 0 line originates 
from the transition to S(l,l) and not to T^(l,l). For the former we 
expect Ae- =E^±J, with 



e 



the exchange energy, whereas the latter would result in Ae^^Ae" 
(in both scenarios, causes an extra fixed offset in both Ae-, but it 
does not contribute to their difference; for transitions to both S(l,l) 
and T^(l,l), we plotted the expected Ae-(B) dependence for our 
DD parameters in the Supplementary Fig. 3a). This interpretation 
is consistent with the increase of Ae^ with larger interdot tunnel 
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Figure 4 | Power-dependence of the photon-assisted tunnelling spectra and spontaneous relaxation, (a) A/gp^ as a function of the double-dot detuning 
eand the microwave amplitude £'for 11 GHz and 6 = 1. OT (middle panel) and 6 = 2. 5T (lower panel). The excitation scheme from Figure Id is employed 
with reference pulse amplitude P8 = 1.5mV. The green (red) dashed lines mark the multi-photon Am = 0 (Am= +1) PAT transitions. The reference 
signal stemming from the pulse is marked by orange arrows. The voltage amplitude E is measured at the end of the coaxial lines at room temperature. 
The uppermost panel displays a linecut measured at 1 T. The red line is a least-squares fit by a sum of four lorentzian peaks to the data, (b) The fitted 
Lorentzian peak area of the Am = +1 line from a is plotted as a function of microwave amplitude E and magnetic field 6. (c) Normalized QPC current 
averaged over t= 5 |is immediately after full mixing at the ST* anticrossing for various magnetic fields at zero microwave power. The spontaneous spin 
relaxation after mixing is measured at a distance Pg to the ST* mixing point (see scheme in Fig. 2c). The dependence on the averaged time interval ris 
displayed in the inset for 6 = 1. 5T and Pg = 2mV together with a least-squares fit (red solid line), (d) The fitted Lorentzian linewidth at half maximum 
(FWHM) of the Am = +1 and the lines Am = 0 from a is plotted as a function of microwave amplitude E'for 6 = 1 T. The error bars in b,d are determined from 
the Lorentzian least-squares fit (see the upper panel of Fig. 4a). 



coupling, hence larger / (Fig. 3a inset; the slopes are not affected). 
In Supplementary Note 3 'The Am = 0 PAT transition, we give more 
arguments for our interpretation by analysing quantitatively the 
intercept Ae^ (B = OT) for various microwave frequencies, as plotted in 
Supplementary Figure S3b. 

So far only single-photon processes were considered, but at 
higher microwave power, also multi-photon lines emerge (Fig. 4a), 
mostly for the S(0,2) -S(l,l) transition (green dashed lines in Fig. 
4a). Like the single-photon S(0,2) -S(l,l) line, their position in 
detuning is B-independent. 

Discussion 

Having shown the power of spin-flip PAT for detailed molecular 
spin spectroscopy, we now discuss the mechanisms responsible for 
this process as confirmed by our simulations. As a first possibility, 
the transitions from S(0,2) to the triplet (1,1) states can take place 
through a virtual process involving S(l,l): the state S(0,2) is coupled 

to S(l,l) by the interdot tunnel coupling, and an (effective) magnetic 

field gradient = (AB^, 0, Ab" ) across the DD couples the spin part 
of all the (1,1) states to each other Here AB has a contribution 
from the effective nuclear field and from the [i magnet. The transi- 
tion matrix element from S(0,2) to T-(l,l) is °^ t^{AB^ / B) , assum- 
ing J. In the following, we use the B-dependence of this process 
as a fingerprint and focus on the red Am = +1 line in Figure 2a,b, as 
we can follow it over the entire magnetic field range. The intensity 
of this line is constant in B, and even if the microwave amplitude E 
is varied, we observe no B-dependence in the area under this peak 
(Fig. 4b). This does not in itself provide evidence that the transition 
rate is magnetic field independent. 



As stated above, the observed PAT lines reflect the steady- state 
change in the charge configuration resulting from stimulated 
photon emission and absorption, and spontaneous relaxation. A 
field-independent steady state could thus be reached from a field 
dependence of relaxation and excitation that cancel each other. We 
therefore verify that the spontaneous relaxation rate is field-inde- 
pendent as well. The measurement of the spontaneous relaxation 
rate is done as follows: Starting from S(0,2), we populate the T+(l,l) 
by 50% through a 200-ns detuning pulse^^ with amplitude Pe in the 
absence of microwaves, and monitor the decay back to S(0,2). The 
relaxation rate can be extracted from the time-averaged lock-in 
signal AJqp(2(t) = (AJq /t) exp(-tr^)dt, where Tis the time spent 
in Pauli blockade between the pulses. AIo = Iqpc(0), which is inde- 
pendent from By is extracted from the fit in the inset of Figure 4c. 
To cover various values of the detuning and the magnetic field, we 
next fix T=5|is and record Algp^ as a function of for three dif- 
ferent magnetic fields. We observe that A/qpc(5 jis) and thus are 
essentially independent of B (Fig. 4c). This holds true for all P^ and 
hence for all T+(l,l)-S(0,2) energy splittings. Note that regardless 
of B, this energy splitting is given by Pe alone. This reflects exactly 
the situation in the PAT experiment, for which the microwave fre- 
quency alone sets the energy splitting and, thus, also the required 
phonon energy for the spontaneous relaxation process. The field- 
independence of relaxation suggests that the coupling mechanism 
is magnetic field independent, and, thus, virtual processes involving 
S(l,l) do not give a strong contribution to the transition rates. 

More recently, two mechanisms were considered that provide 
a direct, 5-independent matrix element between S(0,2) and the 
(1,1) triplet states: (i) the hyperfine contact Hamiltonian is of the 
form S .8{r - rj)Ij - S . Thus, nuclear spins, Ij , in the barrier regions, 

5 
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Table 1 1 Microwave transitions in the (1,1) spin manifold. 



Transition Q (MHz) y(|is-^) 



|70^|S(1,1)> 5.3 1.9 

\r)^\T^) 1.5 0.16 

\r)^\T-) 1.0 0.04 

|S(T1)>^|70> 27 2.0 

|S(1,1)>^|7-> 18 1.9 

|PH|7-> 5.1 0.20 



The Rabi frequencies Q and dephasing rates y, for transitions between different pairs of states 
of the (1,1) spin manifold, are calculated with the model explained in Methods (detuning 
e=-40|ieV, 6g^t=1Tand lOjieV microwave amplitude). The dephasing rates include only 
phonon decay and co-tunnel effects. Further spin dephasing due to, for example, hyperfine 
coupling is neglected. 



where the left- dot and right- dot orbitals overlap, can flip-flop with 
the electron spin, S, simultaneously with charge tunnelling^^ (ii) 
The SO Hamiltonian is of the form p^^yS^^ and can directly couple 
states that differ in both orbital and spin^*^'^^ (When the orbital part 
of the initial and final state are the same, the SO Hamiltonian does 
not provide a direct matrix element, and the transition rate becomes 
B-dependent.^'^^"^^) The ratio of the SO -mediated rate and the hyper- 
fine mediated rate can be estimated as [( A ViV / 4A)(a/Aso)]^»whichis 
a few thousand in the experiment (see Supplementary Note 4 'Spin 
flip tunnelling mechanism' for the derivation). Here A is the single- 
dot level spacing, iVthe number of nuclei in contact with one dot, A 
the hyperfine coupling strength, a the interdot distance and X^^ the 
spin-orbit length. We, therefore, believe that SO interaction is the 
dominant spin-flip mechanism for the observed PAT transitions. The 
presence of a magnetic-field-independent matrix element between 
S(0,2) and the (1,1) triplets is confirmed by the observation that the 
intensity of the reference signal shows no field dependence. 

Finally, we extract from Figure 4a the fitted linewidth as a func- 
tion of driving power (Fig. 4d). For small £, we find a width ^ 1 GHz, 
similar to that observed earlier for spin- conserving PAT processes^. 
For stronger driving, both the Am = 0 and the Am = + 1 lines are 
further broadened, up to £~ 1.5 mV. If these lines were power broad- 
ened, their width would imply transition rates in excess of 1 GHz. 
However, we do not believe that this is the case, because, in meas- 
urements with short microwave bursts, the populations saturated 
only on a long (10-|is) timescale (data not shown). Presumably, 
charge or gate-voltage noise is responsible for the broadening 
instead. 

We have shown that in our DD system, all (1,1) spin states have 
at least weakly allowed electric-dipole transitions to S(0,2). In mate- 
rials with high SO interaction like In As, the effect of the non-spin 
conserving PAT will be even stronger. In materials with weak SO 
interaction, a strong gradient magnetic field can be used to facili- 
tate spin-flip tunnelling transitions. In all cases, all (1,1) to S(0,2) 
transitions can be driven. This opens up the possibility of realizing 
full quantum control of the two spins in the (1,1) manifold through 
off- resonant (microwave) -stimulated Raman transitions through 
the excited S(0,2) state, at small negative detuning and in the pres- 
ence of a finite ABU. Furthermore, we find from our analysis that 
the vicinity in energy of the S(0,2) state also allows direct (single- 
frequency) transitions to be induced between any two -spin states 
within the (1,1) subspace. Direct and stimulated Raman transitions 
are governed by the same matrix elements (see Methods). 

The calculated Rabi frequencies Q,^ for the direct transitions are 
shown in Table 1. For each pair of states, the dephasing rate Jij from 
phonon decay and co-tunnelling effects is calculated as well. In the 
absence of spin dephasing, rotations will be limited by the achiev- 
able ^ijljip with the overall amplitude damping in time t going as 
tYijl2, leading to an infidelity of a TT-pulse of order 7r'X./(2Q..). Con- 
sidering dephasing due to the random, quasi- static nuclear field 
with characteristic dephasing time T2 -10-30 ns, we find for our 



device parameters that 7/7^1/'?2- Therefore, phonon mediated 
decay does not increase dephasing beyond the hyperfine contribu- 
tion. Importantly, a quasi- static dephasing contribution (here from 
the hyperfine field) affects Rabi oscillations less strongly than rap- 
idly fluctuating fields (here from phonons); it permits Rabi oscilla- 
tions to be observed with periods much longer than T2 ^^'^^ 

Controllability of the ( 1 , 1 ) manifold can be best achieved by using 
only the three couplings between the triplet manifold and the sin- 
glet, corresponding to the largest Rabi frequencies. With envelope 
amplitudes of 10 |aeV, just within the perturbative limits of our 
theory, Rabi frequencies of order 10 MHz are possible, similar to the 
frequencies obtained in similar quantum dots for single-spin rota- 
tions. Thus, we find that the complete (1,1) subspace is controllable 
using direct transitions near the S(0,2) crossing, with many rota- 
tions possible^'^^'^^ enabling microwave -induced entangling gates. 

In addition, our observations suggest new measurement tech- 
niques that do not rely on Pauli spin blockade^^. An example is a 
measurement that distinguishes parallel from anti-parallel spins 
while acting non- destructively on the S(l,l) - T^(l,l) subspace, 
which can be achieved by coupling resonantly the T^(l,l) and 
T"(l,l) states to S(0,2), followed by charge readout. This constitutes 
a partial Bell measurement and leads to a new method for produc- 
ing and purifying entangled spin states that enables universal meas- 
urement-based quantum computation^^"^^ 

Methods 

Sample fabrication. 30-nm thick TiAu gates are fabricated on a 90-nm deep 
(AIq 3,Gao 7)As/GaAs two-dimensional electron gas (2DEG) by means of e-beam 
lithography. The double-dot axis is aligned along the [110] GaAs crystal direction 
(z-direction), which is parallel to the external magnetic field direction B (Fig. 
la). The 2DEG is Si 5-doped (40 nm away from the hetero-interface), exhibits an 
electron density of 2.05xl0^^cm"^ and a mobility of 2.06xl0^cm^ per Vs at 1 K in 
the dark. The grounded, 275 nm thick, 2-\im wide and 10-|im long Cojl magnet 
is evaporated on top of an 80-nm thick dielectric layer, aligned along B (magnetic 
easy axis) and placed -400 nm away from the closest dot centre. We calculate^^ 
that, at the double-dot position, the |i magnet adds ^^-l 10 mT to B^ and generates 
a magnetic field gradient of ABU = 5 niT per 50 nm and a transverse gradient of 
A5^~ - 6 mT per 50 nm at saturation {B>2 T). 

Measurement. The sample is mounted in an Oxford KelvinOx 300 dilution re- 
frigerator at 30 mK. Left and right side gate voltages, V^l ^i^d V^, are set by low-pass 
filtered dc lines and -60 dB attenuated coaxial lines combined with bias-tees with 
a cutoff frequency of 30 Hz. The pre-amplified current through the quantum-point 
contact is read out by a lock- in amplifier locked to the 880 Hz on-off modulation of 
the microwaves. The bias across the double dot is set to 0|iV. Voltage pulses to the 
left and right side gates are generated with a Sony Textronix AWG520. The micro- 
waves are generated with a HP83650A and combined with the pulses to the right 
side gate. Microwave bursts and detuning pulses are synchronized to ensure that 
the microwave excitation is switched off during the detuning pulses that generate 
the reference signal (Fig. Id). 

The gates constitute an open-ended termination of the transmission line and, 
therefore, the microwaves predominantly generate an AC electric field giving rise 
to purely electric-dipole PAT transitions. The AC magnetic field contribution from 
the displacement current is negligible at the dot position. The highest contribution 
is the displacement current from the right side gate to the 2DEG. The separation 
between them is 90 nm and the shortest distance between the right dot and the 
displacement current is more than 100 nm. The relevant area of the displacement 
current is the area of the right side gate closest to the dot, which is roughly 40 nm 
by 200 nm. At an AC voltage of 1 mV and at the highest frequency of 20 GHz, the 
displacement current generates an AC magnetic field of less than 2 nT at the right 
dot, corresponding to a Rabi frequency of 1 1 Hz. This is three orders of magnitude 
slower than the spin relaxation time (Fig. 4) and therefore negligible (see also the 
reasoning in ref. 9). 

Simulation. The. Hamiltonian describing the two-spin system near the 
(1,1) - (0,2) transition is taken to be a five-state system, with four (1,1) spin 
states and a (0,2) spin singlet^^ in the presence of an external magnetic 
field B and a magnetic field gradient AB , which includes both the quasi- 
static nuclear field and the field from the |i magnet. This is given by 
H = g/lBPuiA ■ (B + AB) + S2-{B- AB))Pu - £ \ S(0,2)XS(0,2) | -FH^, were 
P^i is the projector onto the (1,1) subspace, e is the detuning due to the 
difference in gate potentials from the left and right gates and the tunnel coupling 
H,= a|S(l,l)XS(0,2)|) + U|r^(l>l)X5(0>2)|) + U|7^-(l>l)X5(0,2)|) + h.cwitht, 
the spin-conserving tunnel coupling and t^^ the spin-orbit coupling set to 5% of t^. 
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To find the signal, we expect theoretically from the experiment, we add a weak, 
rapidly oscillating term to the Hamiltonian: £-^£q + Qcos(vt). We diagonalize H 
with Q = 0, then make a rotating frame transformation in which levels are grouped 
into bands n (defined by a projector P„) where the states in a band n are much 
closer in energy than /zv, while the energy difference between states in band n and 
« + 1 are within 2/3rds of hv. Each band rotates at a rate nv, and we can then make 
a rotating wave approximation, keeping terms due to d that couple adjacent bands, 
that is, our perturbation in the rotating frame and rotating wave approximation is 
V = Q 2^ PJ (0, 2)SX(0, 2)S I P^+i + h.c. 

Next, we add dissipation and dephasing by including relaxation due to 
coupling of the electron charge to piezoelectric phonons in a two -orbital 
(Heitler-London-like) model. We thereby neglect deformation phonons as 
the energy scales examined in the experiment (7-22 GHz) are much smaller 
than the characteristic frequency scale of a phonon on the length scale of the 
dot c^;,//jof~60-120 GHz. To determine the coupling, we take as an ansatz for 
the electronic wavefunctions the Fock-Darwin states, given by Gaussians, and 
calculate the coupling after orthogonalizing the states with the perturbation 
^ph = ^e,kf(K ) / P^e,k y^ 'Pe -^l-k^' where/(/cJ~l for the energy 
scales we are working with. We then use Fermi's golden rule to calculate excitation 
and relaxation from thermal and spontaneous emission of phonons. Finally, we 
numerically solve the superoperator for the steady state and compare the expecta- 
tion value of |S(0,2))(S(0,2)| with and without the excitation Q., mimicking the 
effect of the lock-in detection. 

Coherent transitions within the (1,1) manifold. We conclude with a discussion of 
the controllability of the (1,1) charge manifold using modulated microwaves on e. 
To evaluate this, we consider the high- magnetic field, large detuning limit, where 
-£,glilB \^tg ^ gjil^ I AB I, where is the characteristic tunnelling energy 
scale. In this limit, we can use a Schrieffer- Wolff transformation for Hq + y into the 
(1,1) subspace, where 



(2) 



in the basis (r-(l,l), |iT>, |Ti>, T+(l,l), S-(0,2)) with b = gii^B^,db = gji^AB^^ . 
Defining the un-normalized coupling vector 

\v) = u\T~) + t |iT> - + u\T^), the tunnel coupling is represented by 

V =\ v){S{0, 2)\ + \ S{0, 2))(v I , with u = tso,y , t = Htso,z - ^5 ) / ^/2 . We seek S 

such that exp(A5)(Ho + Ay)exp( - A5) is block diagonal (has no terms coupling (1,1) 
to (1,1)) to order A^. Here A is a perturbation theory parameter that we will take 
to be unity after the analysis. S is given by solution of P( [S,Ho] + y)Q = 0, where P 
projects onto the (1,1) charge space and Q projects onto the (0,2) configuration. It 
takes the form S = |5> <S(0,2)| - |S(0,2)^ <5|, with 

|5> = ^|r-> + ^|iT>-^|n> + ^|r-> 

s + b £-db £ + db s-b 

Finally, for the (1,1) manifold, we get the effective Hamiltonian 
Hgff =HQ+^{\B}{v\ + \v}{B\to order A^ of Hq +1/2[S,V]. 

Within this framework, we can now evaluate the effects to order A^ expected 
from a time-dependent perturbation of the form = Q{t)cos{vt)\S{0,2)){S{0,2)\. 
We get coupling between (1,1) and (0,2) at order A, and coupling between (1,1) 
states at order A^. The charge-changing coupling becomes (taking A^l): 

e^^ I S(0, 2)XS(0, 2) I e'^^ =\ S{0, 2)XS(0, 2) | +X[S, \ S(0, 2)XS(0, 2)\] + 
^[S,[S,|S(0,2)XS(0,2)|]] 



fb 
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0 
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db 
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-b 
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.0 
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0 





=\B){S{0,2)\ + \S{0,2)}{B\ + \BXB\ + 
{1-{B\B))\S{0,2)XS{0,2)\. 



(3) 



From this expression, we see that microwave transitions are possible through 
time- dependent driving of e, if |5> has a non-zero matrix element for the desired 
state. We also see immediately that Raman coupling through |S(0,2)> has the same 
matrix elements as a direct coupling between (1,1) states through the |5>(5| term. 

Thus, full controllability of the (1,1) manifold through microwaves can be 
arrived at by two methods. One is to use stimulated Raman transitions through the 
|S(0,2)> state, the other is to use direct transitions at order A^. Both yield similar 
coupling strengths so that the stimulated Raman process does not beat the direct 
process in this case. However, we remark that in the absence of a magnetic field 
gradient, controllability of the |TnI>, |slT> manifold, without coupling to the |r-> 



states, can not be achieved this way, as there is zero coupling through microwaves 
within this theory to the pure | T°> state. However, with the magnetic- field gradients 
in the experiment, controllability is recovered with the constraint that operation 
pulses cannot be faster than db/h. 
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